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ABSTRACT 


2015 RAS 


We present the identification of potential members of nearby Galactic globular 
clusters using radial velocities from the RAdial Velocity Experiment Data Release 4 
(RAVE-DR4) survey database. Our identihcations are based on three globular clus¬ 
ters - NGC 3201, NGG 5139 (w Gen) and NGG 362 - all of which are shown to have 
|RV| > 100 km s“^. The high radial velocity of cluster members compared to the 
bulk of surrounding disc stars enables us to identify members using their measured 
radial velocities, supplemented by proper motion information and location relative to 
the tidal radius of each cluster. The identification of globular cluster stars in RAVE 
DR4 data offers a unique opportunity to test the precision and accuracy of the stel¬ 
lar parameters determined with the currently available Stellar Parameter Pipelines 
(SPPs) used in the survey, as globular clusters are ideal testbeds for the validation of 
stellar atmospheric parameters, abundances, distances and ages. For both NGC 3201 
and uj Gen, there is compelling evidence for numerous members (> 10) in the RAVE 
database; in the case of NGC 362 the evidence is more ambiguous, and there may be 
significant foreground and/or background contamination in our kinematically-selected 
sample. A comparison of the RAVE-derived stellar parameters and abundances with 
published values for each cluster and with BASTl isochrones for ages and metallicities 
from the literature reveals overall good agreement, with the exception of the apparent 
underestimation of surface gravities for giants, in particular for the most metal-poor 
stars. Moreover, if the selected members are part of the main body of eac h cluster our 
results would also suggest that the distances from Binney et al. (2014), where only 
isochrones more metal-rich than -0.9 dex were used, are typically underestimated by 
^ 40 % with respect to t he published distances for the clusters, while the distances 
from Zwitter et al. (2010) show stars ranging from 1 to ~ 6.5 kpc - with indications of 
a trend toward higher distances at lower metallicities - for the three clusters analysed 
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1 INTRODUCTION 

In the era of massive stellar spectroscopic surveys, auto¬ 
mated Stellar Parameter Pipelines (SPPs) and their valida¬ 
tion are crucial for the scientific exploitation both of exist- 


certainty in the distances to globular clusters is ~ 6%, which 


ing Galactic surveys, such as SEGUE (Yanny et al. 20091 


and RAVE (Steinmetz et al. 20061, and those in progress, 
such as Gaia, which will measure spectra for ~ 150 million 


stars (|de Bruijne||2012|), Gaia-ESO ( 

Gilmore et al.||2012|), 

APOGEE ([Allende Prieto et al.[[2008 

), and GAL AH ( 

Free- 

man & Bland-Hawthorn|2008 Zucker et al.|2013| [Anguiano 

et al. 2014 De Silva et al. 

2015 

), where observations of a 


million stars are planned. However, there are limited oppor¬ 
tunities for checking the outputs of these automated SSPs 
against more traditional analyses in the literature, aside 
from dedicated observations of reference or calibration stars. 

Galactic globular cluster (hereafter GGG) members of¬ 
fer a unique opportunity to validate the precision and ac¬ 
curacy of fundamental stellar atmospheric parameters ob¬ 
tained using currently available SPPs ( |Lee et S(] |2008[ 


Smolinski et al. 20111. The GGC population in the halo 


of the Milky Way covers a wide range of metal abundances, 
essentially independent of radius from the Galactic Gentre, 
spanning approximately -0.5 dex to -2.2 dex. In addition 
these objects have an age spread of 2 - 3 Gyrs (Marfn-Eranch 


|et al.|200^|VandenBerg et al.|20f^ , being mostly older than 


10 Gyr (Harris 19961. However, the traditional paradigm 


treating GGGs as single stellar populations has largely fallen 
by the wayside in recent years. Multiple generations of stars 
have been detected from photometry and spectroscopy in 
a number of GGGs ( Piotto et al.|[2007 Milone et al.|[Ml0 


[Gratton et al.|2012[ ). For a number of massive star clusters, 
like u> Gen, several distinct episodes of star formation have 


been discovered (Johnson & Pilachowski 20101. Very mas¬ 


sive star clusters are considered as possible cores or nuclei 


of stripped dwarf galaxies (Bekki & Freeman 20031. High- 


resolution spectroscopic studies of individual stars in GGGs 
have revealed that some of these objects have a substantial 


star-to-star metallicity scatter. Marino et al. (20111 reported 


a range in metallicity from -2.0 dex to -1.6 dex in a data-set 
of 35 red giants in M22. [Carretta et al.| ( |MTo| ), using high- 
resolution spectra of 76 red giants, found that the bulk of 
stars peak at [Fe/H] ~ -1.6, with a long tail extending to 
higher metallicities, in the globular cluster M54. Very precise 
abundance determinations for several members of NGC 3201 
with high-resolution spectroscopy show a possible metallic¬ 
ity spread of 0.12 dex is present in the cluster ([Munoz e^ 


al. 2013). However, except for variations in their light ele¬ 


ment abundances (Garretta et al. 2010), most GGGs seem 


to be mono-metallicity objects; that is, they have roughly 
uniform iron abundances. A real spread in metallicity seems 
to be rare. 

The period of time over which chemically distinct multi¬ 
ple generations of stars in GGGs are believed to have formed 
- ~ 10® to ~ 10® years - ( Conroy Sz Sperg^|2011 Gratton 


et al. 20121 is still one to two orders of magnitude shorter 


than the history of star formation in the Galactic disc, and 
these stars formed within a relatively small volume (r ~ tens 
of pc). Hence any identified GGC members in the RAVE 
survey database would serve as excellent test subjects for 
validating the estimated distances and ages in the RAVE 
survey using stellar atmospheric parameters; the typical un¬ 


leads to a 13% uncertainty in the absolute ages (Sarajedini 
VandenBerg et aL||2013). In this paper, we report on 


2009 


the identification of members of nearby globular clusters in 
the RAVE catalogue, and use these identifications, in con¬ 
junction with the properties of these clusters published in 
the literature, to test the basic stellar properties obtained for 
these stars from the RAVE survey. (An independent search 
of RAVE data for GGC members, with the goal of using 
new extra-tidal stars as tracers of the clusters’ disruption 
and possible accretion origins, has been carried out by jKun-j 
der et al.|2014 


This paper is organised as follows. In Section 2 we de¬ 
scribe the RAVE survey. Cluster membership selection is 
described in Section 3. In Section 4 we use the likely cluster 
members to test the stellar parameters in the RAVE survey. 
We present our conclusions in Section 5. 


2 THE RAVE SURVEY 


The RAdial Velocity Experiment (RAVE; [Steinmetz et al.] 
20061 is a spectroscopic survey which used the Six De¬ 


gree Field (6dF) multi-object spectrograph on the 1.2m UK 
Schmidt Telescope of the Australian Astronomical Observa¬ 
tory (AAO) at Siding Spring Observatory (SSO). Stars were 
initially drawn from the pilot survey input catalogue based 
on Tycho-2 ( [Hpg et al.|2000[ ) and SuperCOSMOS ( [Hambl}^ 
et al.[ 20011 before the main input catalogue based on DE- 
NIS ( [Epchtein et al.|1997 l, in the magnitude range 9 < I < 
13. The RAVE survey provides, via the fourth data release 
(DR4), the radial velocities and stellar atmospheric param¬ 
eters for 483,849 objects, derived using medium-resolution 
spectra (R = 7,500) in the Ca H triplet region (8410 - 8795 
A). RAVE data are complemented with proper motions from 
Tycho-2, UGAG2, UGAC3, PPMX, PPMXL and SPM4, in 
addition to photometric data from the major optical and 
infrared catalogs, Tycho-2, USNO-B, DENIS, 2MASS and 
APASS (Steinmetz et al.[[2006[ Zwitter et al.[[2008 Siebert 


et al. 2011 Kordopatis et al. 2013). Boeche et al. (2011) 


presented the RAVE chemical catalogue. It contains chemi¬ 
cal abundances for the elements Mg, Al, Si, Ga, Ti, Fe, and 
Ni, with a mean error of ~0.2 dex, as judged from accuracy 
tests performed on both synthetic and real spectra. |Breddel^ 
et al. (2010) developed a method for estimating distances 


from RAVE spectroscopic data, stellar models and (J-Ks) 
photometry from archival sources to derive absolute magni¬ 
tudes. Zwitter et al. ( 2010| determined new distances with a 
method assuming that the star undergoes a standard stellar 
evolution and that its spectrum shows no peculiarities. [Bur- 


nett et al. (2011) applied the Bayesian scheme of Burnett 


& Binney (2010) to the DR3 data. Recently, Binney et al. 


(2014) utilised the Burnett & Binney (2010) technique on 


the DR4 parameters with the addition of H-band photome¬ 
try and a determination of the extinction to estimate stellar 
distances. 


3 THE SAMPLE 

We identified potential stellar members for three GGGs us¬ 
ing the RAVE catalogue. The clusters are NGC 3201, NGC 
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Figure 1. Galactic coordinates plot of all the stars observed in the RAVE survey. The large red dots indicate the positions of globular 
clusters for which we identified potential members in the survey database. 


5139 {uj Cen) and NGC 362. We note that we also have 
identified potential stars from numerous other clusters - in¬ 
cluding NGC 2298, NGC 2808, NGC 4833, NGC 5897, NGC 
6496, NGC 6541, NGC 1904 (M79) and NGC 1851 - in the 
RAVE database, but, as the systemic velocities of these clus¬ 
ters are not sufficiently separated from the velocity distri¬ 
bution of Galactic foreground stars, we cannot determine 
likely membership without taking abundances into account, 
thereby introducing a metallicity bias into our selection cri¬ 
teria. Hence we restrict our analysis in this study to the 
aforementioned three kinematically distinct GGCs. 

Table 1 summarises the main properties of each of the 
clusters included in this study. Positions and radial velocities 
come from the compilation of the Harris (19961 catalogue 
(2010 edition) while the absolute proper motions of globular 
clusters come from Casetti-Dinescu et al. (2007 2010 20131 
and references therein. Fig. 1 shows the positions of all the 
stars observed in the RAVE survey in Galactic coordinates. 
The red points indicate the positions of the GGCs found in 
this study. 



0 10 20 30 40 50 60 

error_RV (km/s) 


Figure 2. RAVE radial velocity measurement errors with respect 
to the S/N of the spectrum. For this study we selected stars with 
S/N > 10 (indicated by a solid line in the plot), where ~ 80% 
have a ajiy ^ 5 km s“^. 


3.1 Cluster membership selection 

We selected probable cluster members based on their mea¬ 
sured radial velocities (RVs) and absolute proper mo¬ 
tions. RAVE RVs are computed by cross-correlation of sky- 
subtracted normalised spectra with an extensive library of 
synthetic spectra. Spectra without sky subtraction are used 


to compute the zero-point correction. The typical RV ac¬ 
curacy for RAVE data is ^ 2 km s“^. For measurements 
with a high signal-to-noise ratio (S/N), the error is only 1.3 


km s , with a negligible zero-point error (Steinmetz et al. 


2006[ [Zwitter et al.|2008| [Siebert et al.|2011f [Kordopatis ^ 


al.|26i3|. 

In Fig. we present the behaviour of the error in the 
RV with the S/N of the spectra. (For more detail on S/N 
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Table 1. Properties from the [Harris ] | |1996[ | catalogue (2010 edition) of GGCs identified in the RAVE data with |RV| > 100 
km s“^. Distance uncertainties are consistent with a ± 0.15 dex change in the distance modulus. 


ID 

RA(hh:mm:ss) 

(J2000) 

DEC(dd:mm:ss) 

(J2000) 

1(°) 

b n 

Vr (km/s) 

Pa (mas/yr) 

PS (mas/yr) 

Dist. 

(kpc) 

NGC 362 

01:03:14.26 

-70:50:55.6 

301.53 

-46.25 

±223.5 ± 0.5 

±4.873 ± 0.514 

-2.727 ± 0.824 

8.6 ± 0.6 

NGC 3201 

10:17:36.82 

-46:24:44.9 

277.23 

±8.64 

±494.0 ± 0.2 

±5.280 ± 0.320 

-0.980 ± 0.330 

4.9 ± 0.3 

NGC 5139 

13:26:47.24 

-47:28:46.5 

309.10 

±14.97 

-1-232.1 ± 0.1 

-5.080 ± 0.350 

-3.570 ± 0.340 

5.2 ± 0.3 


Table 2. Structural parameters of Galactic Globular Clus¬ 
ters identified in the RAVE data from the [Harris 1 l |1996| l 
catalogue (2010 edition). 


ID 

Central ay 
(km/s) 

Tc 

(arcmin) 

rt 

(arcmin) 

NGC3201 

5.0 ± 0.2 

1.30 

25.34 

NGC5139 

16.8 ± 0.3 

2.37 

48.38 

NGC362 

6.4 ± 0.3 

0.18 

10.35 


brightest giants of a cluster. Fig. shows the RVDF for a 6 
square-degree area centred on each cluster, in which we have 
identified potential members in the RAVE data. We made 
a broad selection between dwarfs and giants using RAVE- 
derived surface gravities and temperatures (which, owing to 
their different luminosity classes, probe different volumes). 
The black line represents stars in the range 3500 K < Tefi 
< 8000 K and log g > 3.5 (cgs) while the red line represents 
stars in the range 3600 K < Teff < 6000 K and log g < 3.5 
(cgs). 


measurements in the RAVE spectra see Zwitter et al.|2008 ) 
In this work we selected stars with S/N > 10 (black line in 
figure); these stars have a typical internal RV error ^ 5 km 
s“^. Note that RAVE spectra with low S/N (i.e. S/N < 10) 
and stars with spectral peculiarities ( Matijevic et al.||2012 l 
typically have lower precision RVs. Fig. shows that the 
typical accuracy of the RV increases with the S/N, with error 
in RV < 2 km s“^ for 90% of the catalogue at S/N ~ 100. 
We also selected stars with a Tonry & Davis ( |1979| ) cross¬ 
correlation coefficient R larger than 5. RAVE provides very 
precise RVs for 80% of the catalogue. For this reason our 
main criterion for membership selection was a comparison 
of the stellar RVs in the area of the sky around each GGC 
to the systemic RV for that GGC from the literature (see 
Table 1). 

In this paper we restricted our study to the three clus¬ 
ters that have |RV| > 100 km s“^, in order to get cleaner 
samples of halo stars relatively uncontaminated by disc 
stars. This purely kinematically derived sample does not re¬ 
quire us to make any cuts in metallicity for membership 
selection, thereby leaving abundances as a free parameter 
for directly testing the reliability of the RAVE SPP. Once 
candidates were selected as explained above we used the 
tidal radius of the clusters and the proper motions from the 
PPMXL catalogue ( Roeser et al.|20l0 | to select the highest 
likelihood members for the final dataset. 


3.1.1 Selection of potential members based on radial 
velocities 

Nearby GGCs that have a large RV (|RV| > 100 km s“^) 
can serve as unbiased abundance calibrators for stellar sur¬ 
veys, as these objects can be identified in the RV Distribu¬ 
tion Function (RVDF hereafter) from the bulk of the Milky 
Way’s stellar disc for the area in the sky surrounding the 
cluster. RAVE is a magnitude-limited survey covering 9 to 
13 in the I band. The three clusters we consider here are a 
few kpc away; in this magnitude range, the candidate mem¬ 
bers of one of these clusters will likely belong to the upper 
parts of the red and asymptotic giant branches, i.e., the 


NGC 3201 has an extremely high radial velocity (RV 
~ -1-494 km s“^), which makes this cluster an ideal target 
to test our selection method. High-velocity halo star studies 
using the RAVE survey have shown that the number of stars 
with a RV larger than ±400 km s“^ is not very large | Smith 


et al.|2007 Piffi et al.|2014 1. In Fig. 3 we identify a group of 


eleven stars with a RV ~ 495 km s“^ in the area of the sky 
around NGC 3201. This group of stars is clearly distinct in 
velocity from the bulk of the stars in this region of sky, and, 
given their similar RV values, they are very likely members 
of NGC 3201. The central velocity dispersion for this cluster 
is only 5 km s“^ (see Table 2), which is greater than the 
mean RV error (< 3 km s“^). 


Candidates for the globular cluster NGC 5139 (a; Gen) 
also appear clearly in the RAVE data. In the RVDF there 
is a peak of stars around 230 km s“^, exactly the RV re¬ 
ported for this cluster (see Table 1). As mentioned above, 
uj Gen is a massive, complex system with multiple star for- 


mation episodes (Joo & Lee 

2013|l. Da Costa & Coleman 

(2008 

1, Wylie-de Boer et al. 

(20101 and Majewski et al. 


(20121 detected field stars that may be associated with uj 
Cen in the nearby Galactic disc. We found a broad RV dis¬ 
tribution around the nominal RV for this cluster, from ~ 
200 to 260 km s“^ (Fig. 3); however, it is important to note 
that this cluster has a signific ant central i nterna l dispersion, 
av = 16.8 km s“^. Recently, Da Costa (20121 measured a 
line-of-sight velocity dispersion in the outer parts of uj Cen 
of ~ 6.5 km s“^. 


The case for identifying members of NGC 362 is less 
clear-cut in the RAVE data. This cluster has an RV ~ ±223 
km s“^. The RVDF in Fig. 3 shows only three potential 
candidates for NGC362 membership. We also find a poten¬ 
tial stellar radial velocity over-density between 100 < RV < 
200 km s“^; some of these stars 
Small Magellanic Cloud (SMC) 

Tg^. 


may be associated with the 
(Da Costa & Hatzidimitriou 
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NGC3201 


NGC5139 




NGC362 



Figure 3. Radial Velocity Distribution Functions (RVDFs) for an area of 6 square degrees centred at the positions of NGC 3201, NGC 
5139 (a; Cen) and NGG 362. The black line represents dwarf stars and the red line giants (see text for an explanation of the selection 
of dwarfs and giants). Note the peak around 490 km s“^ related with NGG 3201, around 230 km s“^ associated with the RV for NGC 
5139 and around 230 km s“^ for NGC 362. The distribution of stars between 120 < RV < 200 km s“^ in the panel for NGC 362 is 
intriguing; they may be associated with the Small Magellanic Cloud (SMC). 


3.1.2 Selection of probable members using tidal radii and 
proper motions 

As a consistency check to verify the validity of our probable 
cluster members selected using RVs only, here we check the 
proper motions of the candidates as well as their position 
within the clnsters’ tidal radii. 

The central potential of the clnster dehnes the density 
distribution within the cluster, where the radius is in units 
of the King radius, ro (KJngjlQMl. At the radius where the 
potential falls to zero, the density also falls to zero. This is 
the so-called tidal radius, rt, and it is the radius at which 
the inward force towards the cluster centre is balanced by 
the outward pull of the Galaxy’s tidal field, so that at r > rt 
stars are not bound to the cluster. The central concentration 
parameter of a King model is defined to be 


c= logic (n/r-o) (1) 

where ro is the radius at which the projected surface 
brightness falls to half its central value. Table 2 lists the core 
radius and the tidal radius for the globular clusters presented 
in this study, as well as the internal velocity dispersion of 
the clusters. 

For the cluster NGC 3201 we selected eleven stars using 
RV alone. Nine of these stars are inside the tidal radius (see 
Fig. 13) with proper motions similar to the nominal value 


for the cluster {nc, cos{5) = 3.65, ns = 0.98). We also find 
four stars with significantly higher proper motion values, 
suggesting that these stars may be foreground objects. 

Table summarises the candidates of NGC3201. The 
two stars outside of the tidal radius are J102025.9-46440610 
and J101405.6-46284110; however both of these stars have 
similar proper motions to the cluster, and hence are still can¬ 
didates for membership. The stars with high proper motions 
(J101725.9-462621, J101738.6-462716, J101751.5-46221010 
and J101716.2-462533) could be catalogued as field stars 
despite being inside the tidal radius. Despite the marked 
differences in proper motion compared to NGC 3201, the 
striking similarity in RV makes these stars interesting tar¬ 
gets. A detailed chemical abundances analysis is necessary 
to understand their relation, if any, with the cluster. All the 
stars have spectra with a S/N > 30, except for two with S/N 
> 20 . 

We found a total of 52 stars to be potential members 
of u! Cen based on their RVs (Tablej^). 31 of those stars are 
inside the tidal radius and the proper motion diagram in 
Fig[^ shows that 45 stars share an apparent proper motion 
close to the value for oj Cen (/Tc cos(5) = -3.43, fj ,5 = -3.57), 
suggesting these stars are probable members of oj Cen. 

10 stars were selected using the RV for the cluster NGC 
362 (Table [ 3 . We find five stars inside the tidal radius (see 
Fig. §. These stars are very close to each other but, un- 
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Table 3. NGC 3201 candidates selected from RAVE data and their parameters. The last column indicates the crowding (see Section 3.2). 


ID 

R.A. 

Deck 

RV 

(km/s) 

pa COs((5) 
(mas/yr) 

(mas/yr) 

^J■5 

(mas/yr) 

(mas/yr) 

J 

S/N 

C. 

J101405.6-462841 

lOh 14’ 05.63” 

-46° 28’ 40.6” 

494.56 ± 0.62 

1.76 

12.00 

3.20 

12.00 

10.01 

35.9 

- 

J101640.5-463221 

lOh 16’ 40.52” 

-46° 32’ 20.5” 

493.94 ± 0.83 

-0.33 

12.00 

12.00 

12.00 

10.14 

34.5 

U 

J101648.9-461807 

lOh 16’ 48.88” 

-46° 18’ 06.9” 

496.44 ± 0.91 

-0.82 

12.00 

7.40 

12.00 

9.16 

33.7 

- 

J101716.2-462533 

lOh 17’ 16.16” 

-46° 25’ 32.8” 

502.19 ± 2.57 

74.46 

12.00 

54.50 

12.00 

9.26 

45.7 

V, u 

J101725.9-462621 

lOh 17’ 25.86” 

-46° 26’ 21.0” 

489.22 ± 1.05 

27.03 

14.10 

47.30 

14.00 

8.69 

23.1 

V, u 

J101731.6-462901 

lOh 17’ 31.59” 

-46° 29’ 01.0” 

500.99 ± 1.67 

-14.91 

12.00 

-11.20 

12.00 

10.67 

22.0 

- 

J101738.6-462716 

lOh 17’ 38.59” 

-46° 27’ 16.0” 

494.65 ± 1.41 

129.70 

32.00 

-267.20 

32.00 

- 

36.5 

V, u 

J101751.5-462210 

lOh 17’ 51.53” 

-46° 22’ 09.7” 

495.87 ± 0.68 

37.01 

13.00 

-6.70 

13.00 

9.47 

46.8 

V, u 

J101752.1-461407 

lOh 17’ 52.04” 

-46° 14’ 06.6” 

483.86 ± 0.61 

-10.78 

12.00 

-15.50 

12.00 

9.71 

50.6 

- 

J101859.1-463438 

lOh 18’ 59.10” 

-46° 34’ 37.6” 

497.94 ± 0.86 

-4.18 

3.00 

0.10 

3.01 

8.44 

43.0 

- 

J 102025.9-464406 

lOh 20’ 25.86” 

-46° 44’ 05.8” 

496.77 ± 0.69 

-2.68 

12.00 

10.10 

12.00 

8.41 

40.8 

- 



276.0 276.5 277.0 277.5 278.0 278.5 

Galactic longitude (°) 


100 
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a -100 
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-200 


-300 


Figure 4. Top panel: The positions (in Galactic coordinates) for 
the stars in the RAVE catalogue in the region of NGC 3201. The 
red dots are potential members of NGC 3201 selected based on 
RV. The black circle is the tidal radius according to the value 
given in ||Harris| ( |1996| l (2010 edition) for this cluster. Bottom 
panel: The proper motion plane (mas/yr) for the selected can¬ 
didates. The red dot indicates the absolute proper motion of the 
cluster ( [Casetti-Dinescu et al.|2007[|2010[|2013| |. 
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fortunately, we do not have proper motion information for 
them. The remaining stars are clearly outside of this radius. 
Fig.i shows that the proper motion of these stars is large 
with respect to the cluster {fj,a cos{S) = -0.79, = -2.73). 


There is one star, J005038.4-732818, with proper motions 
close to zero and a similar RV to the NGC 362, that might 
be associated with the cluster despite being outside of the 
tidal radius. 

In summary, using the RVs and proper motions for stars 
in the RAVE catalogue and cluster tidal radii from the lit¬ 
erature we were able to identify - without any abundance 
selection bias - likely candidate members for several nearby 
galactic globular clusters. 


3.2 Crowding 


The 6dF spectrograph used for the RAVE survey placed 
wide (6.7 arcsec) hbres on the sky; the size of these h- 
bres could present problems because of the potential for 
overlapping background spectra of other stars, especially in 
the dense stellar regions around globular clusters. Fibres in 
which the composite spectra of more than one star overlap 
should not be used for the calibration of stellar parameters, 
as the accuracy of the resulting analyses can be at best un¬ 
certain. To identify targets with possible crowding problems 
we have done two different tests. Stars classified as prob¬ 
lematic with respect to general crowding or the presence 
of close neighbours via a visual inspection of DSS/2MASS 
finder charts are marked with “v” in Tables ii and[^ We 
also used UCAC4 ( Zacharias et al.|2013 l to mark suspicious 
cases; namely, we checked if there were any UCAC4 stars 
within 9 arcsec, or any relatively bright stars at larger sep¬ 
arations. The stars for which this is true are marked with 
“u” on the tables. In many cases, the crowding problem also 
led to uncertain (large) proper motions in the PPMXL cat¬ 
alogue, especially in the case of NGC 362. 


4 STELLAR PARAMETERS IN THE RAVE 
SURVEY 

The RAVE survey uses the wavelength region AA8410 - 8795 
A for the determination of the main stellar parameters in 
the atmosphere of the observed stars. In this region we find 
the Ca II triplet, iron and a-element lines; these features are 
prominent spectral lines even for spectra with low signal-to- 
noise and for metal-poor stars, making this region useful for 
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Table 4. NGC 5139 {lo Cen) candidates selected from RAVE data and their parameters. The last column indicates the crowding. 


ID 



R.A. 



Decl. 


RV 

(km/s) 

[let cos((5) 
(mas/yr) 

(mas/yr) 

(mas/yr) i 

(mas/yr) 

J 

S/N 

I 

c. 

J131313.9-460352 

13h 

13’ 

13.88” 

-46° 

03’ 

52.0” 

180.59 

± 

0.84 

11.79 

11.00 

6.00 

11.00 

10.33 

45.1 


- 

J131340.4-484714 

13h 

13’ 

40.43” 

-48° 

47’ 

13.8” 

220.25 

± 

1.07 

-1.58 

12.00 

-4.50 

12.00 

10.63 

33.0 


- 

J131511.1-455458 

13h 

15’ 

11.06” 

-45° 

54’ 

57.5” 

165.54 

± 

1.88 

-2.63 

11.00 

47.60 

11.00 

10.65 

14.8 

V, 

u 

J131548.1-443935 

13h 

15’ 

48.07” 

-44° 

39’ 

34.7” 

314.76 

± 

0.59 

7.09 

1.00 

-10.40 

1.00 

9.50 

57.3 


- 

J131602.3-480507 

13h 

16’ 

02.28” 

-48° 

05’ 

06.6” 

257.20 

± 

0.95 

8.59 

9.00 

11.80 

9.00 

9.66 

21.1 


- 

J131613.1-452004 

13h 

16’ 

13.13” 

-45° 

20’ 

03.7” 

238.33 

± 

1.53 

-2.73 

2.00 

-6.80 

2.00 

10.72 

35.8 


- 

J131729.6-462521 

13h 

17’ 

29.62” 

-46° 

25’ 

21.3” 

196.70 

± 

0.93 

1.22 

11.00 

1.10 

11.00 

10.63 

42.5 


- 

J 132045.6-445053 

13h 

20’ 

45.61” 

-44° 

50’ 

52.9” 

300.91 

± 

1.15 

4.02 

1.00 

-5.20 

1.00 

10.83 

44.3 


- 

J132209.4-481432 

13h 

22’ 

09.41” 

-48° 

14’ 

31.7” 

175.63 

± 

0.62 

3.45 

9.00 

-0.30 

9.00 

8.79 

54.4 


- 

J132430.7-472427 

13h 

24’ 

30.74” 

-47° 

24’ 

26.5” 

229.33 

± 

0.75 

13.44 

11.00 

1.70 

11.00 

8.92 

41.5 


- 

J 132446.8-472449 

13h 

24’ 

46.76” 

-47° 

24’ 

48.6” 

234.78 

± 

1.10 

2.87 

11.00 

4.10 

11.00 

10.51 

42.0 


u 

J132517.5-472427 

13h 

25’ 

17.53” 

-47° 

24’ 

26.6” 

228.57 

± 

0.57 

1.82 

3.00 

-5.10 

3.00 

9.54 

57.1 


- 

J132521.3-473654 

13h 

25’ 

21.32” 

-47° 

36’ 

54.0” 

239.65 

± 

0.40 

7.32 

3.00 

-5.50 

3.00 

8.66 

91.9 


- 

J132521.8-452320 

13h 

25’ 

21.79” 

-45° 

23’ 

19.8” 

162.52 

± 

0.62 

-1.86 

12.00 

11.20 

12.00 

10.45 

36.2 


- 

J132545.2-473238 

13h 

25’ 

45.15” 

-47° 

32’ 

38.3” 

217.80 

± 

0.68 

2.46 

3.00 

-7.10 

3.00 

9.39 

49.9 

V, 

u 

J132551.2-472702 

13h 

25’ 

51.20” 

-47° 

27’ 

01.8” 

222.19 

± 

0.49 

9.47 

3.00 

-5.30 

3.00 

9.18 

47.8 

V, 

u 

J132552.0-473016 

13h 

25’ 

52.00” 

-47° 

30’ 

16.3” 

213.89 

± 

0.89 

-58.49 

3.00 

-5.50 

3.00 

8.81 

52.9 

V, 

u 

J132558.7-473610 

13h 

25’ 

58.72” 

-47° 

36’ 

09.8” 

235.57 

± 

2.20 

7.01 

3.00 

-4.00 

3.00 

9.11 

34.7 


- 

J132601.7-474034 

13h 

26’ 

01.74” 

-47° 

40’ 

33.6” 

253.80 

± 

0.65 

5.79 

3.00 

-0.20 

3.00 

9.27 

66.4 

V, 

u 

J132609.1-472720 

13h 

26’ 

09.06” 

-47° 

27’ 

19.5” 

236.99 

± 

0.55 

5.59 

2.00 

-9.40 

2.00 

9.66 

51.2 

V, 

u 

J132614.6-472123 

13h 

26’ 

14.56” 

-47° 

21’ 

22.9” 

244.65 

± 

0.60 

11.38 

3.00 

-8.70 

3.00 

8.67 

44.2 


V 

J132623.7-474243 

13h 

26’ 

23.66” 

-47° 

42’ 

42.5” 

233.22 

± 

1.00 

-1.33 

8.00 

1.70 

8.00 

9.19 

30.1 


- 

J132629.6-473701 

13h 

26’ 

29.62” 

-47° 

37’ 

01.4” 

231.27 

± 

1.13 

2.82 

3.00 

-9.90 

3.00 

9.35 

59.0 


- 

J 132639.0-474359 

13h 

26’ 

38.96” 

-47° 

43’ 

58.5” 

226.20 

± 

1.22 

0.16 

3.00 

-3.20 

3.00 

10.05 

41.2 


- 

J132639.3-472035 

13h 

26’ 

39.27” 

-47° 

20’ 

34.8” 

216.40 

± 

0.64 

-68.91 

12.00 

27.50 

12.00 

9.24 

66.4 

V, 

u 

J132646.2-471415 

13h 

26’ 

46.17” 

-47° 

14’ 

15.2” 

223.66 

± 

0.66 

4.27 

2.00 

-4.80 

2.00 

9.29 

63.2 


- 

J 132654.3-474605 

13h 

26’ 

54.34” 

-47° 

46’ 

05.1” 

240.49 

± 

0.64 

2.56 

3.00 

-11.10 

3.00 

10.01 

80.3 


- 

J 132704.4-443003 

13h 

27’ 

04.43” 

-44° 

30’ 

03.1” 

328.80 

± 

1.32 

-5.99 

9.00 

26.40 

9.00 

11.08 

24.9 


V 

J132709.6-472052 

13h 

27’ 

09.60” 

-47° 

20’ 

51.5” 

236.96 

± 

0.40 

10.53 

2.00 

-6.80 

2.00 

9.05 

60.6 

V, 

u 

J132710.5-473701 

13h 

27’ 

10.54” 

-47° 

37’ 

00.5” 

248.80 

± 

1.11 

-89.77 

18 

-86.90 

18.00 

9.08 

44.4 


V 

J132726.0-473060 

13h 

27’ 

25.97” 

-47° 

30’ 

59.9” 

261.82 

± 

1.19 

5.45 

3.00 

-8.30 

3 .00 

9.37 

59.0 

V, 

u 

J132753.7-472442 

13h 

27’ 

53.72” 

-47° 

24’ 

42.1” 

234.62 

± 

1.86 

-0.09 

13.00 

-15.50 

13.00 

10.61 

38.3 

V, 

u 

J132754.7-471932 

13h 

27’ 

54.71” 

-47° 

19’ 

32.2” 

248.04 

± 

2.23 

6.47 

2.00 

-3.00 

2.00 

9.08 

24.5 


- 

J132757.3-473638 

13h 

27’ 

57.31” 

-47° 

36’ 

38.1” 

252.55 

± 

0.97 

70.18 

15.00 

18.30 

15.00 

10.58 

42.2 

V, 

u 

J 132800.8-473247 

13h 

28’ 

00.77” 

-47° 

32’ 

47.0” 

248.72 

± 

0.64 

-1.00 

3.00 

-9.60 

3.00 

10.34 

67.9 

V, 

u 

J 132804.8-474504 

13h 

28’ 

04.83” 

-47° 

45’ 

04.1” 

220.21 

± 

0.75 

4.86 

2.00 

-9.60 

2.00 

8.65 

34.3 


- 

J132813.6-472424 

13h 

28’ 

13.57” 

-47° 

24’ 

23.5” 

231.31 

± 

1.04 

3.74 

9.00 

-1.90 

9.00 

9.41 

30.0 


- 

J132815.0-473739 

13h 

28’ 

15.04” 

-47° 

37’ 

39.4” 

237.78 

± 

1.88 

28.64 

11.00 

-8.00 

11.00 

9.94 

20.0 


- 

J132816.9-472956 

13h 

28’ 

16.90” 

-47° 

29’ 

56.0” 

245.12 

± 

0.56 

8.18 

3.00 

-8.30 

3.00 

9.02 

55.4 


- 

J 132833.8-473206 

13h 

28’ 

33.81” 

-47° 

32’ 

05.6” 

239.87 

± 

0.79 

-1.10 

3.00 

-9.90 

3.00 

8.84 

64.7 


- 

J 132839.9-472633 

13h 

28’ 

39.93” 

-47° 

26’ 

32.9” 

228.40 

± 

2.89 

5.51 

3.00 

-9.20 

3.00 

8.82 

23.9 


V 

J132918.9-471924 

13h 

29’ 

18.86” 

-47° 

19’ 

23.9” 

233.41 

± 

0.64 

1.86 

9.00 

-0.90 

9.00 

9.01 

60.0 


- 

J 132936.8-500005 

13h 

29’ 

36.76” 

-50° 

00’ 

05.4” 

199.04 

± 

0.83 

-0.58 

2.00 

3.60 

2.00 

8.03 

44.3 


- 

J133106.0-483312 

13h 

31’ 

06.04” 

-48° 

33’ 

11.9” 

205.74 

± 

0.92 

1.79 

11.00 

5.10 

11.00 

10.25 

41.1 


- 

J133257.2-492045 

13h 

32’ 

57.15” 

-49° 

20’ 

45.2” 

182.31 

± 

0.67 

2.60 

10.00 

7.70 

10.00 

9.71 

55.8 


V 

J133328.7-441903 

13h 

33’ 

28.71” 

-44° 

19’ 

02.8” 

223.66 

± 

0.65 

8.97 

2.00 

-7.10 

2.00 

7.73 

74.8 


- 

J133430.0-474615 

13h 

34’ 

29.96” 

-47° 

46’ 

14.8” 

161.65 

± 

1.49 

5.74 

3.00 

-11.80 

3.00 

8.89 

21.9 


- 

J133536.3-450533 

13h 

35’ 

36.28” 

-45° 

05’ 

33.0” 

168.34 

± 

0.79 

-7.91 

2.00 

-2.00 

2.00 

8.41 

58.7 


- 

J 133609.4-440854 

13h 

36’ 

09.39” 

-44° 

08’ 

53.6” 

181.41 

± 

1.30 

-7.99 

1.00 

-19.20 

1.00 

9.80 

29.3 


- 

J133628.5-461932 

13h 

36’ 

28.54” 

-46° 

19’ 

32.1” 

162.48 

± 

0.97 

3.42 

12.00 

-1.60 

12.00 

9.45 

30.7 


- 

J133939.4-490014 

13h 

39’ 

39.41” 

-49° 

00’ 

13.5” 

236.59 

± 

0.56 

0.79 

2.00 

-35.00 

2.00 

9.66 

61.5 


- 

J134508.6-492832 

13h 

45’ 

08.64” 

-49° 

28’ 

32.1” 

217.44 

± 

0.68 

-3.59 

12.00 

-11.10 

12.00 

8.52 

68.8 


- 


metallicity estimations over a broad range of stellar proper¬ 
ties (Cenarro et al. [2001 Kordopatis et al.|2011|. 


The methodology for determining stellar parameters in 
the su rvey has undergone several revisions. [Zwitter et al.| 
120081 used a penalised 'V method employing an extensive 


grid of synthetic spectra calculated from the latest version 


of Kurucz stellar atmosphere models for the first and sec¬ 
ond RAVE data releases. From comparison with external 
data sets, [Zwitter et al.| ( |2008| ) estimated errors in stellar 
parameters for a RAVE spectrum with an average signal-to- 
noise ratio (S/N ~ 40) to be 400 K in temperature, 0.5 
dex in gravity, and 0.2 dex in metallicity. jSiebert et al.j 
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Table 5. NGC 362 candidates selected from RAVE data and their parameters. The last column indicates the crowding. 


ID 

R.A. 



Deck 

RV 

(km/s) 

pa COs((5) 
(mas/yr) 

(mas/yr) 

Ts 

(mas/yr) 

(mas/yr) 

J 

S/N 

C. 

J004905.3-733108 

OOh 49’ 05.27” 

CO 

o 

31’ 

07.9” 

229.27 ± 0.78 

-6.59 

2.00 

-8.00 

2.0 

9.20 

31.2 

- 

J004217.1-740615 

OOh 42’ 17.14” 

-74° 

06’ 

15.3” 

266.93 ± 0.79 

5.01 

1.00 

-15.10 

1.0 

9.40 

35.5 

- 

J005038.4-732818 

OOh 50’ 38.39” 

-73° 

28’ 

18.3” 

221.15 ± 7.19 

0.29 

1.00 

-1.90 

1.00 

11.07 

36.0 

- 

J010313.6-705037 

Olh 03’ 13.62” 

O 

O 

50’ 

36.8” 

220.97 ± 1.09 

- 

- 

- 

- 

- 

58.6 

V, u 

J010314.7-705115 

Olh 03’ 14.67” 

o 

O 

51’ 

15.3” 

225.40 ± 1.48 

- 

- 

- 

- 

- 

26.4 

- 

J010314.7-705059 

Olh 03’ 14.74” 

o 

o 

50’ 

58.9” 

225.01 ± 0.55 

- 

- 

- 

- 

- 

106.1 

V, u 

J010315.1-705032 

Olh 03’ 15.10” 

o 

O 

50’ 

32.3” 

233.88 ± 1.10 

- 

- 

- 

- 

- 

40.6 

V, u 

J010319.0-705051 

Olh 03’ 19.03” 

o 

O 

50’ 

51.4” 

222.78 ± 1.01 

- 

- 

- 

- 

- 

51.4 

V, u 

J010335.7-705052 

Olh 03’ 35.71” 

o 

O 

50’ 

52.0” 

228.24 ± 0.74 

7.94 

2.2 

58.6 

2.2 

9.73 

38.6 

V, u 

J011655.9-690607 

Olh 16’ 55.87” 

-69° 

06’ 

07.3” 

236.51 ± 1.80 

11.80 

10.00 

-7.00 

10.00 

9.92 

18.9 

- 



Galactic longitude (°) 



IIq, cos( 5) (mas/yr) 


Figure 5. Top panel: The positions (in Galactic coordinates) for 
the stars in the RAVE catalogue in the region of to Gen. The 
red dots are potential members of lj Cen selected based on RV. 
The black circle is the tidal radius according to the value given in 
|Harris| | |1996[ | (2010 edition) for this cluster. We identify several 
stars outside of the tidal radius as candidate cluster members. 
Bottom panel: The proper motion plane (mas/yr) for the selected 
candidates. The red dot indicates the absolute proper motion of 
the cluster ( [Casetti-Dinescu et al.|2007||2010[[^013[ |. 


(20111, for the third data release, used new synthetic spec¬ 
tra for intermediate metallicities that were added in order 
to provide a more realistic spacing toward the densest re¬ 
gion of the observed parameter space, and thereby remove 


biases toward low metallicity. They also used a new con¬ 
tinuum normalisation which signicantly reduced the corre¬ 
lation between metallicity and S/N, masked bad pixels and 
improved the radial velocity zero-point. Finally, |Kordopatis| 


et al. (20131 computed the stellar atmospheric parameters 


using a new pipeline, based on the algorithms of MATISSE 
( Recio-Blanco et al. |200^ [Kordopatis et al.|2011 1 and DE¬ 
GAS (DEcision tree alGorithm for Astrophysics) for the 
fourth RAVE data release (DR4). Spectral degeneracies and 
2MASS photometric information are also taken into consid¬ 
eration. In this study we use the DR4 stellar parameters, as 
derived in Kordopatis et al. (20131. Boeche et al. (20111 pre¬ 
sented elemental abundances derived from RAVE spectra for 
the elements Mg, Al, Si, Ga, Ti, Fe, and Ni, through a pro¬ 
cessing pipeline in which the curve of growth of individual 
lines is obtained from a library of absorption line equivalent 
widths to construct a model spectrum, that is then matched 
to the observed spectrum via a x^-minimisation technique. 

In this section we use a subset of the candidate cluster 
members identified in the previous section, namely those 
stars which match the radial velocities and fall within the 
tidal radii of the three globular clusters discussed above - 
NGG 3201, u! Cen and NGC 362 - to test the validity of the 
stellar parameters derived from RAVE spectra. Since most 


GGCs appear to have ages > 10 Gyr (e.g. Dotter et al. 
|2010[ [VandenBerg et al.|2013P we decided to use isochrones 
for the ages determined for these clusters in the literature to 
test Teft, [m/H] and log g determinations. Although to Cen 
may contain younger stellar populations (e.g., |Villanova et| 
"^pOOTt , it is very difficult to derive the ages of single gi¬ 


ant stars to within a few Gyr given that temperatures along 
computed red giant branches (RGBs) are very sensitive to 
many aspects of stellar physics (notably convection) and be¬ 
cause the location of the RGB is much less dependent on age 
than on metal abundances (e.g., [VandenBerg et al.||2012| ). 
We also explore the abundances for the elements Mg, Al, Si, 
Ca, Ti, Fe, and Ni for the potential members of the clusters. 


4.1 The case of NGC 3201 


NGC 3201 shows very peculiar kinematic characteristics. 
The cluster has an extreme radial velocity and a highly ret¬ 
rograde orbit ( Casetti-Dinescu et al. [2007 1. Thus kinemati¬ 
cally NGC 3201 appears likely to be of extragalactic origin; 
however Munoz et al. (20131 claim that its chemical evolu- 
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Figure 6. Top panel: The positions (in Galactic coordinates) for 
the stars in the RAVE catalogue in the region of NGC 362. The 
red dots are potential members of NGG 362 selected based on RV. 
The black circle is the tidal radius according to the value given 
in [Harris ] (2010 edition) for this cluster. We identified a 

number of stars outside of the tidal radius as candidate members 
of the cluster. Bottom panel: The proper motion plane (mas/yr) 
for the selected candidates. The red dot indicates the absolute 
proper motion of the cluster (jCasetti-Dinescu et al.||2007| |2010| 
|2013| |. 


tion was similar to most other, presumably “native”, GGCs. 
NGC 3201 is a low mass halo cluster and the existence of 
star-to-star metallicity variations remains controversial. In 
agreement with the findings of|Da Costa et aI~|(|1981[),|Car- 


retta et al. (20091 found no significant variations of [Fe/H] 


in NGC 3201 in their analysis of high-resolution spectra of 
hundreds of stars, and [Munoz et al.| | |2013| | similarly found 
no evidence for any intrinsic Fe abundance spread, except 
for one star in their sample. On the other h and, jGonzalezj 
jfc Wallerstein j ( |1998| ) and jSimmerer et al.~| ( [2013 1 identi¬ 
fied a spread in [Fe/H] in cluster stars at least as large as 
0.4 dex, even though [Covey et al. (20031 were not able to 
confirm the presence of a significant spread in metallicity 
within the cluster greater than about 0.3 dex. Iron abun¬ 
dance variations have been found in the most massive glob¬ 
ular clusters (e.g., M22, M54, uj Gen), suggesting multiple 
star formation episodes and metal enrichment via Type la 
supernova events. Simmerer et al. ] (2013l concluded that a 
real [Fe/H] spread, if it did exist, would support the idea that 
NGC 3201 was initially far more massive, formed outside the 


Milky Way, and was subsequently captured. Recently, a pos¬ 
sible solution to the controversy surrounding the metallicity 
spread in NGC 3201 has been proposed by jMucciarelli et al.j 
(20151, who demonstrated that the metal-poor component 
claimed by Simmerer et al. (20131 is composed of asymp¬ 
totic giant branch stars that could be affected by non-local 
thermodynamical equilibrium effects driven by iron overion¬ 
isation, and therefore concluded that there is no evidence of 
intrinsic iron spread. 


4.1.1 NGC 3201 abundances, temperatures and gravities 


RAVE DR4 metallicities were calibrated using [Fe/H] from 
the literature and dedicated observations of calibration 
stars. Since the metallicity measurement is dominated by 
the Ca H lines, we have the overall metallicity [M/H]Dfl 4 ~ 
[Fe/H] -b a small correction from a-elements. In this work 
we use [M/H]_d_r 4 = [Fe/H], see Fig. 8 in Kordopatis et al. 
( 2013||. For this c luster we adopt <[Fe/H]> = -1.5 ± 0.1 dex 
(Carretta et al. 1 (200^, which agrees well with other deter- 


20131 


minations (e.g. Gonzalez fc Wallerstein ||1998 Munoz et al. 


From the highest likelihood members of NGC 3201 - 
i.e., those with RVs and proper motions consistent with 
membership and which also fall within the tidal radius - we 
find two groups of stars (see Fig.j^ black dots). The largest 
group has a mean [Fe/H] of-1.80 ±0.11 dex while the second 
group contains only two stars and has a <[Fe/H]> ~ -1.30 
± 0.10 dex. The RAVE SPP is able to identify these stars 
as metal-poor and the mean [Fe/H] value for this cluster 
from the RAVE pipeline is -1.55 ± 0.10, in excellent agree¬ 
ment with high-resolution spectroscopy as indicated above. 
The first group has an observed scatter of 0.05 dex, a lower 
value than the nominal uncertainty of the metallicity mea¬ 
surement, suggesting that there is no evidence for an in¬ 
trinsic scatter in the metallicity in terms of RAVE errors. 
The second group exhibits an observed scatter of 0.07 dex, 
and again there is no evidence for an intrinsic scatter, al¬ 
though there are only two stars in the group. Thus we have 
found two clearly different metallicity groups that are po¬ 
tentially members of NGC 3201 based solely on their kine¬ 
matics and proper motions. However these two groups could 
be just an artefact of selection effects on our small sample 
size, as other authors with larger data-sets have not reported 
such a bimodal distribution. Combining the two groups we 
find an observed scatter of 0.24 dex for NGC 3201, larger 
than the expected metallicity uncertainties. RAVE metal¬ 
licities suggest that an intrinsic scatter is observed in this 
cluster with a star-to-star metallicity variation from -1.25 
to -1.83 dex, however the metallicity of the two metal-rich 
stars should be verified before we draw conclusions on the 
intrinsic scatter in this cluster. This spread in metallicity 


is in good agreement with Gonzalez & Wallerstein (19981 


and Simmerer et al. (20131; however Simmerer et al. (20131 


found their stars ranged between -1.80 < [Fe/H] < -1.40, and 
[Gonzalez fc Wallerstein (19981 found a spread over -1.65 < 
[Fe/H] < -1.15. As mentioned above, some authors did not 
find a spread in metallicity in this cluster. Moreover, the 
published colour-magnitude diagrams for the RGB of NGC 
3201 do not appear to show enough scatter to accommodate 
a range in [Fe/H] between -1.3 and -1.8 dex, in particular 
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Figure 7. Left panel: Radial velocity-abundance plot for RV candidate members which fall within the tidal radius of NGC 3201. The 
dashed lines indicate the nominal RV and [Fe/H] for the cluster reported in the literature. The blue dots are stars with large proper 
motions with respect to the measured proper motion for the cluster. Right panel: The metallicity distribution function for all RV candidate 
members which fall within the tidal radius (in red) and for those which also have proper motion measurements consistent with that of 
the cluster (black line). 


the CMDs corrected for differential reddening (von Braun 
fc Mateo|2001| . 


There are four stars with proper motions significantly 
different from the nominal proper motion for the cluster. 
These stars are shown in blue in Fig. Two stars lie in the 
most metal-poor group, in good agreement with the over¬ 
all metallicity; one star is in the metal-rich group but it 
is slightly more metal-rich; and one star has an intermedi¬ 
ate metallicity. Could these stars with high proper motions 
in fact be halo stars and not members of NGC 3201? This 
scenario seems rather unlikely. There are not many stars 
with extreme heliocentric RVs (> 400 km s”'^) in the Galac¬ 
tic halo (e.g., Smith et al. 2007| Piffl et ^|2014 although 
it should be noted that the RVs in these two studies are 
Galactocentric and not heliocentric). For NGC 3201, RVcaj 
~ -1-275 km s“^; this velocity is still high, but both papers 
identihed objects in this velocity regime in the RAVE cat¬ 
alogue. These stars also have a metallicity range in good 
agreement with the abundances reported for this cluster, 
and they are at projected distances close to the cluster. 
However, typical halo stars can have metallicities in this 
range. Other potential explanations for the proper motion 
discrepancies are that these stars could be members of NGC 
3201 but significantly closer to us than the main body of the 
cluster itself, or that the proper motions for these objects 
are simply imprecise or erroneous as these objects are in 
crowded fields. 


Harris (19961 reported that the stars in NGC 3201 are 


around 10.5 Gyr old. Simmerer et al. (20131 found that 14 


of the observed stars but they concluded that a younger 
age for the metal-poor stars would improve the isochrone 
ht. [Potter et al.| ( |2010| ) derived an age of 12.0 ± 0.75 Gyr 
from a deep HST CMD of NGC 3201, assuming that the 
stars in the cluster are a-enhanced. Recently, [Munoz et al.j 
( 2013[ ) reported an age of 11.4 Gyr, taking into account the 
(G-l-N-l-0) abundance in the isochrone age. 


In addition, [Munoz et al.[ ( [2013| ) found, using only Mg, 
Si, Ga and Ti, a mean a-element abundance for stars in 
NGG 3201 of -1-0.30 ± 0.06. This value is in good agreement 
with the [a/Fe] values estimated by [Gonzalez fc Wallerstein[ 


(19981 and Garretta et al. (20091. In order to test the 


Gyr isochrones with [a/Fe] = 0.0 fit the stellar parameters 


effective stellar temperatures and the surface gravities de¬ 
rived from the RAVE spectra we compared them to BASTI 
isochrones ( Cassisi et al.|200'6 1 . We selected isochrones with 
an age = 11.5 Gyr and [Fe/H] = -1.31 and -1.84, the range 
in metallicity we found for this cluster. We also made use 
of a-enhanced isochrones, in particular, [a/Fe] = -1-0.2 and 
-1-0.4. Fig. 1^ shows the Colour-Magnitude Diagram (CMD) 
using 2MASS J and H bands. From the 2MASS catalogue 
( Cutri et al.|2003 1 we selected stars inside a circle with a ra¬ 
dius of 0.1 degrees centred on NGC 3201. The CMD shows 
the RGB, AGB and the HB of NGG 3201. We find that 
the selected stars from the RAVE catalogue lie in the upper 
part of the RGB indicated by red dots in Fig. the blue 
dots are stars with high proper motion with respect to the 
nominal value of the cluster as indicated above. On the right 
panel in Fig.[^we have the TeU - log g diagram for the stellar 
members of NGC 3201 together with the BASTI isochrones. 
Generally there is a good fit within the errors for the given 
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Figure 8. Left panel: 2MASS J, J-H Colour-Magnitude Diagram (CMD) centred on NGC 3201, with a radius of 0.1 degrees. The large 
dots are the stars selected as being probable members of the cluster based on RV and location within the tidal radius; stars with large 
proper motions relative to the cluster are in blue. Note that most of the stars lie on the RGB of NGC 3201. Right panel: Temperature- 
surface gravity diagram for the stars selected as probable members of NGC 3201, overplotted with BASTI isochrones for an age of 11.5 
Gyr, [Fe/H] = -1.3 and -1.8, and [o/Fe] = -1-0.2 and -1-0.4 dex. The blue dots are stars with large proper motions. 


isochrones. However there is a group of stars for which, in 
the metallicity range reported here, a younger isochrone (~ 
10.5 Gyr) would provide a better fit. If we accept the age of 
NGC 3201 as 11.4 Gyr and [a/Fe] ~ 0.3 dex (Munoz et al. 


2013), we can conclude from this exercise that the gravities 


are generally underestimated with respect to the isochrones. 

Note that there is one star with log g = 3.3, clearly 
outside the main group of stars and the isochrones. This 
star is the most metal-rich ([Fe/H] = -1.16) and has a high 
proper motion. The star is marked as potentially affected by 
crowding problems, and a close look at the spectrum shows 
that the “star” is either two separate stars in the same fibre 
or a spectroscopic binary (SB2). This could easily lead to 
an erroneous measurement of log g. Interestingly, however, 
its estimated metallicity is similar to the other two high 
metallicity stars, neither of which appear to be impacted by 
crowding problems or show composite spectra. 

In the Boeche et al.| ( |2011[ ) chemical catalogue we find 
five candidate stars with [a/Fe] measurements (see Table]^. 
We exclude the one star with [a/Fe] = -1-0.64, as this is the 
object mentioned above, with a high surface gravity and a 
possible composite spectrum. For the other four stars, the 
measured <[a/Fe]> = -1-0.22 is in good agreement with high 
resolution studies, and we find (T[a/_Fe] = 0.18 dex. 


4-1-2 NGC 3201 stellar distances 


To date there have been four studies (Breddels et al.|2010 


Zwitter et al.|20T0 Burnett et al.|[20TT Binney et al.|20r^ 

which address the challenge of calculating distances for the 


RAVE stars. These have primarily used atmospheric stellar 
parameters derived from the spectra, photometric colours 
of the stars and stellar evolutionary tracks, i.e., the method 
of spectrophotometric parallaxes. In this work we test the 


distances derived in Zwitter et al. (2010) and Binney et al. 
(|2014p using the stellar parameters from RAVE DR4 (|Ko 


rdopatis et al. 2013). Zwitter et al. (2010) assumed that 


stars follow standard stellar evolution tracks, as reflected by 
theoretical isochrones. They also assumed that interstellar 
reddening is negligible because the vast majority of RAVE 
stars lie at high Galactic latitudes (] b ] > 20°). With these 
assumptions, Zwitter et al. (2010) determined the probabil¬ 


ity distribution function (PDF) for the absolute magnitude. 
The authors concluded that their derived distances of both 
dwarfs and giants match the astrometric distances of Hip- 
parcos stars ( van Leeuwen|20OT ) to within ~ 21% using the 
RAVE DR3 stellar parameters (Zwitter et al. 2008). Another 


approach, with some assumptions via prior functions, is the 


Bayesian framework developed in Burnett et al. (2011) and 


Binney et al. (2014). Prior functions play a key role in these 


new distances and they reflect the state of our knowledge of 
the Galaxy, with different prior probabilities based on mod¬ 
els of the density of the three components of the Galaxy (thin 
disc, thick disc and halo). Binney et al. (2014) included the 
effects of interstellar dust by applying a prior that reflects 
increasing extinction with distance and higher extinction to¬ 


wards the Galactic plane using the Schlegel maps (Schlegel 
et al.|1998 ) in the prior. The distance determination of Bin 


ney et al. (2014) made use of stellar isochrones which only 
went down to ~ -0.9 dex in metallicity, thus neglecting more 


© 2015 RAS, MNRAS 000, [iH^ 






































































12 B. Anguiano et al. 


Table 6. NGC 3201 candidate members selected from RAVE data and their stellar parameters 


ID 

Te// (K) 


[Fe/H] 

log g (cgs) 

[a/Fe] 

J-H 

J101405.6-462841 

4510 

± 

96 

-1.83 

± 

0.11 

1.17 

± 

0.31 

-1-0.19 

0.72 

J101640.5-463221 

4619 

± 

96 

-1.83 

± 

0.11 

1.05 

± 

0.31 

-1-0.31 

0.64 

J101648.9-461807 

4438 

± 

96 

-1.79 

± 

0.11 

0.70 

± 

0.31 

- 

0.74 

J101716.2-462533 

4797 

± 

96 

-1.16 

± 

0.11 

3.34 

± 

0.31 

-1-0.64 

0.91 

J101725.9-462621 

4248 

± 

96 

-1.50 

± 

0.11 

0.51 

± 

0.31 

- 

0.85 

J101731.6-462901 

4272 

± 

96 

-1.35 

± 

0.11 

0.38 

± 

0.31 

- 

0.82 

J101738.6-462716 

4404 

± 

96 

-1.72 

± 

0.11 

0.16 

± 

0.31 

- 

- 

J101751.5-462210 

4342 

± 

75 

-1.76 

± 

0.09 

0.75 

± 

0.20 

- 

0.79 

J101752.1-461407 

4344 

± 

75 

-1.82 

± 

0.09 

0.72 

± 

0.20 

-1-0.40 

0.73 

J101859.1-463438 

4000 

± 

75 

-1.25 

± 

0.09 

1.00 

± 

0.20 

- 

0.88 

J102025.9-464406 

4066 

± 

96 

-1.69 

± 

0.11 

0.62 

± 

0.31 

-0.03 

0.84 


metal-poor stars. They concluded that the expectation of 
parallax may be the most reliable distance indicator and 
found a good agreement between the expectation values of 
the parallaxes and the values measured by Hipparcos for 
the very few stars in common, especially in the case of hot 
dwarfs. 

In Fig. we show the estimated stellar distances for RV 
and tidal-radius-selected members of NGC 3201 together 
with their measured [Fe/H]. The top plot shows the dis¬ 
tances from |Binney et ^||2014[) and t he bottom shows the 
distances calculated by|Zwitter et al. (20101; Binney et al. 


son of the results of the two distance determination methods 


(20141 estimated the distances for six of the stars that we 


selected as likely members of the cluster, while in |Zwitter| 
et al. (20101 we find estimated distances for 11 such ob¬ 


jects. The blue dots in the hgure are stars with high proper 
motions with respect to the nominal value for the cluster. 
For this study we adopted the true distance modulus from 


(Harris 1996 see also Layden & Sarajedini 20031, (m-M)o 
= 13.45, i.e., 4.9 kpc (red vertical line in Fig. 


et al. (2014) distances for the six members have a small 


range from 2.7 to 3.7 kpc with a mean value of 3.1 kpc. If 
these objects are current members of the cluster and not a 
group of stars disrupted from NGC 3201, our results would 
suggest that the distances for the giant stars we identify 
as likely members of NGC 3201 are underestimated by ~ 
2 kpc. Note that in Fig. we found that the gravities for 
these stars appear to be slightly underestimated with re¬ 
spect to the selected isochrones using the latest estimations 
of age, [Fe/H] and [a/Fe] for this cluster. Underestimated 
surface gravities could potentially affect the distance deter¬ 
mination of the RAVE giants, although from a simplistic 
perspective we would expect lower surface gravities at the 
same temperature to yield higher intrinsic luminosities, and 
hence distances which are too large, rather than too small. 
Zwitter et al. ( 2010|) calculated the distances for 11 


members of the cluster (bottom panel in Fig. [^. The blue 
dots in the hgure are stars with high proper motion values 
with respect to the cluster. The distances found by |Zwit-| 
ter et al. (2010) for these stars range from 1.3 to 6.5 kpc. 


They found a group of stars around 6 kpc, but two of those 
have a high proper motion. The reddening to NGC 3201 is 


rather larger, E(B-V)=0.24 mag (Harris 1996), which may 
contribute to the distance spread, as Zwitter et al. (2010) did 


reveals that Binney et al. (2014) found their sample of stars 


to cover a small range in distances while the stars of |Zwitter| 
|et al.| ( |2010[ ) have a signihcant scatter in distance. 

Recently, a new distance calibration has been applied to 
these stars by P. McMillan, using the methodology of |Binne^ 
|et al.| ( |2014[ ) but incorporating isochrones with metallicities 
extending to lower than [Fe/H] = -0.9. The results of this 
work are consistent with the published distance for this clus¬ 
ter, with a mean distance of 4.2 ± 0.8 kpc for the members 
identihed in this study. 


4.2 The case of NGC 5139 {uj Centauri) 

The putative globular cluster uj Gen is a complex stellar 
system. If it is a globular cluster, it is the most massive 
such cluster in the Galaxy, in which clear evidence of mul¬ 


tiple stellar populations has been detected (Lee et al. 1999 
[Smith et al.|[2000j [Johnson fc Pilachowski [2010[ ). On the 


other hand, there is also evidence that uj Gen could be the 


stripped core of a dwarf elliptical galaxy (Bekki & Freeman 
2003 Bekki fc Norris|2006 ). The very bound retrograde or¬ 


bit supports the idea that the cluster entered the Galaxy as 
part of a more massive system whose orbit decayed through 
dynamical friction (Freeman & Bland-Hawthorn 2002 ). This 


cluster exhibits large star-to-star metallicity variations (~ 
1.4 dex). Several studies have found [Fe/H] ranges from ~ 
-2.1 to ~ -0.7 dex using high-resolution spectroscopy for in¬ 
dividual red giants in the cluster; moreover distinct peaks in 
the iron abundance distribution have been detected, suggest¬ 
ing different star formation episodes (Sollima et al. 2005 Vil 


lanova et al.|2007 Marino et al.|2010 ). Estimates of the time 
periods spanned by these different star formation episodes 


vary widely, however. Villanova et al. (2007 2014) found 


a spread in age of 5 Gyr between the youngest and oldest 
members of the cluster, ranging from 8 to 13 Gyr. In con- 
Hilker et al. (2004); Stanford et al. (2006[) found an 


trast. 


age spread on the order of 2 - 3 Gyr, and [Sollima et ah] 


(2005) reported a small or negligible age dispersion. It is 


important to note that Marino et al. (2012) demonstrated a 


not take reddening into account in their analysis. A compari- 


signihcant variation in the G-l-N-l-0 content among uj Gen’s 
stellar populations, which could easily have an impact on 
these age estimates. 

We selected our highest probability uj Cen members 
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NGC 3201 



distance (kpc) Binney et al. 2014 



distance (kpc) Zwitter et al. 2010 

Figure 9. Distance vs. [Fe/H] for the RV and tidal-radius-selected candidate members of NGC 3201. In the top panel the distances 
are those derived in [Binney et al.] l |2014[ |, in the bottom panel the distances are from [Zwitter et al.j ( |2010l l. Blue dots are stars with a 
large proper motion while triangles indicate stars affected by crowding. The vertical red line indicates the distance of the cluster given 
in [Harris I l|1996[l (2010 edition). 


from the RAVE catalogue using a combination of RVs, 
proper motions and location within the tidal radius. In the 
next section we explore the main stellar parameters for these 
objects. 


4-2.1 Abundances, temperatures and gravities 

We find a large spread in [Fe/H] (~ 2.0 dex) for the se¬ 
lected candidate members of the cluster (see Fig 101. The 


abundances range from approximately solar values ([Fe/H] 
~ 0.0) to metal-poor abundances ([Fe/H] ~ -2.2). The blue 
dots in Fig |10| represent stars with RVs similar to the clus¬ 
ter and which lie inside the tidal radius, but which have 
large proper motions with respect to w Cen. The distribu¬ 
tion of [Fe/H] shows three peaks, at [Fe/H] ~ -0.6, -1.3 and, 
the largest, at -1.8 (Fig |10[ right panel). The large range of 
metallicities observed in the RAVE targets is in very good 
agreement with high resolution spectroscopic studies as indi¬ 
cated above. The RAVE metallicity distribution function we 
obtain is also consistent with several star formation episodes 
in the history of the cluster. 


Johnson & Pilachowski (20101, using 855 red giant 


stars, found that the a elements in oj Cen are generally en¬ 
hanced by ~ -1-0.3 dex and exhibit a metallicity-dependent 
morphology. We find a large spread in [a/Fe] for a given 
[Fe/H] (see Fig|TI| and Table 0 using the values from the 
RAVE chemical abundance catalogue ( Boeche et al.|[?011 |. 
The measured [a/Fe] ranges from -0.23 to -1-0.32 dex with a 
mean error of ^0.2 dex ([Boeche et al.|2011 |. Pancino et al. 
(20021 and Villanova et al.| ( |2007[ ) also found a large spread 
in a-elements for a given [Fe/H]. 



Figure 11. [Fe/H] - [o/Fe] diagram for stars selected as likely 
members of ta Cen, plotting the abundances from the RAVE 
chemical catalogue ( [Boeche et al.|2011[ |. We found a large spread 
in [o/Fe] ranging from -0.23 to -1-0.32 dex, with a measurement 
error of ~ 0.2 dex. 


An age spread from 8 to 13 Gyr between different mem¬ 
bers of the cluster was found by Villanova et al. (20071. In 
Fig. we show a Tefi - log g plot of the high probability 
u) Cen members found in the RAVE survey, with the data 
colour-coded by RAVE [Fe/H]: green dots represent stars be¬ 
tween 0.0 < [Fe/H] < -1.0 (note that there is only one star 
at approximately solar metallicity); the red dots are stars 
with -1.0 < [Fe/H] < -1.6; and the blue dots are stars with 
-1.6 < [Fe/H] < -2.2. We overplotted BASTI isochrones with 
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Figure 10. Left panel: Radial velocity-abundance plot for RV candidate members of NGC 5139 (u) Cen). The dashed line indicates the 
nominal RV for the cluster reported in the literature. The blue dots are stars with large proper motions with respect to the measured 
proper motion for the cluster. Right panel: Metallicity distribution function for all RV candidate members which fall within the tidal 
radius (in red) and for those which also have proper motion measurements consistent with that of the cluster (black line). 



Figure 12. 2MASS J, J-H Colour-Magnitude Diagram (CMD) 
centred on lo Cen, with a radius of 0.1 degrees. The big dots are 
the stars selected for being likely members of the cluster based on 
RV and location within the tidal radius; stars with large proper 
motions with respect to the cluster are shown in blue. Most of 
the stars selected using RV and location clearly lie on or near the 
prominent RCB of ui Cen. 


[Fe/H] = [-0.6, -1.5, -1.8], and [a/Fe] = [-1-0.0, -1-0.2], for ages 


of 8.0 Gyr (top panels) and 13 Gyr (bottom panels). (Fig 10 


appears to show three [Fe/H] peaks in our data, at -0.6, -1.3 
and -1.8 dex, hence onr selection of these isochrones.) There 
is a reasonable match (within the errors) between the log 
g, Tefi and [Fe/H] valnes derived from the RAVE spectra 
and the selected isochrones, at least for the more metal-rich 
stars. However, an isochrone younger than 8.0 Gyr would be 
a better fit, especially for the most metal-poor stars (red and 
blue dots). |Villanova et al.| ( |2007[ ) found that the youngest 
stars in the cluster are around 8.0 Gyr. If this age limit is 
correct, and the most metal-poor stars in lo Gen are not 
somehow younger than their more metal-rich counterparts, 
the isochrones would suggest that the spectroscopic gravities 
for the metal-poor stars are underestimated. Note also that 
for the more metal-rich the underestimation in the spectro¬ 
scopic gravities is less evident. 


4.2.2 RAVE stellar parameters and the \ Johnson & \ 

Pilachowski (2010) study 


We identihed 21 stars in common between the oj Cen can¬ 
didates found in RAVE and the high-resolution, high S/N 
spectra of 855 to Cen RGB members obtained and analysed 
for elemental abundances by Johnson fc Pilachowski] (20101. 
The sample includes nearly all RGB stars brighter than V 
= 13.5 and spans uj Gen’s full metallicity range. 

Effective temperatures (Tes) were determined via the 


empirical (V - K) color-temperature relation from Alonso et 
al. 119991 using the recommended value of E(B - V) = 0.12 


(Harris! 19961. Surface gravity estimates were obtained using 
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Table 7. NGC 5139 candidates selected from RAVE data and their stellar parameters 


ID 

Te// (K) 

[Fe/H] 

log g (cgs) 

[o/Fe] 

J-H 

J131340.4-484714 

4516± 96 

-1.36± 0.11 

0.83± 0.31 

- 

0.74 

J131602.3-480507 

4250± 86 

-0.50± 0.10 

2.00± 0.20 

- 

0.95 

J131613.1-452004 

5201± 103 

-0.84± 0.10 

3.20± 0.23 

- 

0.74 

J132430.7-472427 

4350± 96 

-1.83± 0.11 

O.lOi 0.31 

-1-0.10 

0.48 

J 132446.8-472449 

4590± 96 

-1.81± 0.11 

0.73± 0.31 

- 

0.80 

J132517.5-472427 

4265± 101 

-2.00± 0.10 

0.51± 0.35 

-0.04 

0.48 

J132521.3-473654 

4000± 75 

-1.75± 0.09 

0.50± 0.20 

- 

0.80 

J132545.2-473238 

4073± 96 

-1.45± 0.11 

0.61± 0.31 

-1-0.05 

0.60 

J132551.2-472702 

4497± 101 

-2.00± 0.10 

0.50± 0.35 

- 

0.70 

J132552.0-473016 

4250± 96 

-1.75± 0.11 

0.50± 0.31 

-1-0.13 

0.80 

J132558.7-473610 

4256± 86 

-0.64± 0.10 

1.30± 0.20 

- 

0.52 

J132601.7-474034 

4497± 96 

-1.85± 0.11 

0.30± 0.31 

-0.02 

0.79 

J132609.1-472720 

4493± 101 

-2.00± 0.10 

0.50± 0.35 

- 

0.90 

J132614.6-472123 

4279± 96 

-1.75± 0.11 

0.61± 0.31 

-0.08 

0.67 

J132623.7-474243 

4251± 96 

-1.20± 0.11 

0.60± 0.31 

-1-0.32 

0.82 

J132629.6-473701 

4298± 96 

-1.44± 0.11 

O.Oli 0.31 

-0.23 

0.64 

J 132639.0-474359 

4499± 86 

-0.75± 0.10 

2.50± 0.20 

- 

0.59 

J 132639.3-472035 

4250± 75 

-1.75± 0.09 

0.50± 0.20 

-1-0.11 

0.65 

J132646.2-471415 

4466± 75 

-1.78± 0.09 

0.93± 0.20 

-1-0.11 

0.71 

J 132654.3-474605 

4575± 76 

-2.15± 0.10 

0.73± 0.27 

-1-0.13 

0.72 

J 132704.4-443003 

4881± 96 

-1.51± 0.11 

1.96± 0.31 

-0.06 

0.76 

J132709.6-472052 

4258± 76 

-2.00± 0.10 

0.51± 0.27 

- 

0.76 

J132710.5-473701 

4000± 96 

-1.29± 0.11 

0.42± 0.31 

-1-0.14 

0.81 

J132726.0-473060 

4500± 96 

-1.75± 0.11 

l.OOi 0.31 

- 

0.79 

J132753.7-472442 

4577± 86 

-0.76± 0.10 

0.44± 0.20 

- 

0.74 

J132754.7-471932 

3891± 69 

-0.42± 0.09 

4.51± 0.15 

- 

0.52 

J132757.3-473638 

4669± 101 

-2.09± 0.10 

1.23± 0.35 

- 

0.84 

J 132800.8-473247 

4587± 76 

-2.16± 0.10 

0.81± 0.27 

- 

0.71 

J 132804.8-474504 

4251± 96 

-1.25± 0.11 

2.00± 0.31 

- 

0.79 

J132813.6-472424 

4000± 96 

-1.50± 0.11 

O.OOi 0.31 

- 

0.83 

J132815.0-473739 

4249± 96 

-1.25± 0.11 

l.OOi 0.31 

- 

0.74 

J132816.9-472956 

4311± 101 

-2.04± 0.10 

0.79± 0.35 

-0.07 

0.82 

J 132833.8-473206 

4250± 75 

-1.75± 0.09 

0.50± 0.20 

-0.10 

0.71 

J 132839.9-472633 

4500± 112 

O.OOi 0.11 

2.00± 0.24 

- 

0.66 

J132918.9-471924 

4250± 75 

-1.75± 0.09 

0.50± 0.20 

- 

0.70 

J 132936.8-500005 

4500± 96 

-1.25± 0.11 

2.00± 0.31 

- 

0.61 

J133106.0-483312 

3801± 96 

-l.OOi 0.11 

l.OOi 0.31 

- 

0.68 

J133328.7-441903 

4249± 75 

-1.49± 0.09 

0.51± 0.20 

- 

0.85 

J133939.4-490014 

4689± 62 

-0.99± 0.08 

1.88± 0.14 

-0.01 

0.52 

J134508.6-492832 

3800± 96 

-1.25± 0.11 

0.50± 0.31 

-1-0.03 

0.58 


the photometric temperatures and absolute bolometric mag¬ 
nitudes (Mboi), assuming stellar masses of 0.8 Mq. Chem¬ 
ical abundances were determined through standard equiva¬ 
lent width (EW) analyses using the LTE line analysis code 
MOOG. We compared the temperature, metallicity and sur¬ 
face gravity derived from the RAVE spectra and those de¬ 


rived in Johnson & Pilachowski (20101 for the 21 stars in 


common (see Fig. 141. While there is a general correlation 
in Tejf, we found a systematic offset, such that RAVE tem¬ 
peratures are slightly hotter with respect to the photometric 


ones used in Johnson & Pilachowski (20101. For metallicities 


we also found a rough correlation between both methods, al¬ 
though there are two clearly discrepant stars, with RAVE 
metallicities ~ 1 dex higher than the high-resolution abun¬ 
dances. Finally, surface gravities in the small range where 
these giants lie also appear to be offset, with RAVE esti¬ 
mates mostly in the range log(g) ~0.0 - 1.0 and Johnson & 


Pilachowski (20101 values spanning 0.5 - 1.5. The standard 


deviations of the differences in Teff, [M/H] and log g are 158 
K, 0.3 dex and 0.9 dex respectively. Isochrones are affected 
by many uncertainties in their underlying physics associ¬ 
ated with e.g., treatment of convection, surface boundary 
conditions, etc. We found discrepancies between the stellar 
parameters derived from the RAVE spectra and the BASTI 
isochrones used in this work (see Fig. and Fig. 131. The 
lower right panel of Fig. |14|shows data from|Johnson fc Pi-| 


lachowski (20101 for lo Cen and Campbell et al. (20131 for 


NGC 6752 together with BASTI isochrones with [Fe/H] = 
-1.8 dex and -1.5 dex, respectively. We found good agree¬ 
ment between the isochrones and the data from these two 
independent studies, in contrast to the offset between RAVE 
pipelines results and the BASTI isochrones noted above. 
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Figure 13. Temperature-surface gravity diagram for the stars selected as likely members of NGC 5139 {u) Cen). In the top panels we 
plot BASTI isochrones covering the range in [Fe/H] from -1.8 to -0.6 for [a/Fe] = -1-0.0 (left) and -1-0.2 (right) at an age of 8 Gyr. In 
the bottom panel we use the same set of isochrones, with the same range of metallicities and alpha enhancements, but at an age of 13 
Gyr. The green dots represent stars with metallicities between 0.0 < [Fe/H] < -1.0, the red dots are stars with -1.0 < [Fe/H] < -1.6, and 
the blue dots are those with -1.6 < [Fe/H] < -2.2 . 


4-2.3 Stellar distances 


Here we adopt the true distance modulus from Harris {19961 
(2010 edition), (m-M)o = 13.58, i.e., 5.2 kpc. Fig 15 shows 


the stellar distances for NGC 5139 members together with 
their measured [Fe/H]; the top panel shows the distances 


from Binney et al. (20141 and the bottom panel the dis¬ 


tances calculated by Zwitter et al. (20101. The blue dots 


in the figure are stars with high measured proper motions 
with respect to the nominal value for the cluster as described 


above. Binney et al. (20141 found a wide range of distances 


for these stars, from 1 to 7 kpc. Most of the stars are be- 
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Temperature Metallicity 




(K) RAVE DR4 [M/H] RAVE DR4 


Surface gravity 



Johnson & Pilachowski 2010 Campbell et al. 2013 




Ten{K) T,«(K) 


Figure 14. 1:1 comparison for 21 stars in common between the RAVE DR4 catalogue and the data-set from 
metallicity and surface gravity, respectively. The lower right panel show data from 


Johnson Sz Pilachowski 


Johnson &: Pilachowski 


(2010 I. We compared Teff 

(20101 for uj Cen and [Campbell et al.| | [2013| | for the cluster NGC 6752. BASTI isochrones are also showed (Fe/H] = -1.8, -15 dex 
respectively). We found a good agreement between these data-sets and the isochrones. 


tween 2 and 4 kpc with a peak at ~ 3 kpc. The most metal- 
poor stars show a small dispersion in distances. |Zwitter~et| 
E] poTo| also found stars ranging from 1 to 7 kpc, but 
they show a fairly clear trend between [Fe/H] and distances, 
with metal-rich stars estimated to be closer than metal-poor 
ones. It is worth noting that, despite the wide spread in mea¬ 


sured distances, Zwitter et al. (20101 found a group of stars 


with a mean distance around 5 kpc and [Fe/H] between - 
1.0 and -2.2 dex. This value is in good agreement with the 
distances for cu Cen found in the literature. When applied 
to these stars, the aforementioned new distance calibration 
using isochrones extending to [Fe/H] < —0.9 gives a mean 
distance of 4.2 ± 1.2 kpc. 


There is also evidence of held stars that may be associ¬ 


ated with uj Cen in the nearby Galactic disc (e.g., Wylie-de 
Boer et al.|2010 1. If the distances from Binney et al. (20141 


are correct, that might suggest that we found a former part 
of the cluster that is closer to us than the cluster itself. How¬ 
ever, if these stars are in fact part of the main cluster the 
result would suggest that the distances from [Binney et al.| 
(20141 are underestimated by ~ 40%. Note that using the 


distances from Zwitter et al. (20101, we hnd that most of 


the stars are between 4 and 6 kpc, in good agreement with 
distances determined for this cluster from the literature. 


4.3 The case of NGC 362 

NGC 362 has a metallicity [Fe/H] = -1.26, according to the 
2010 update of the Harris] ( [l996| ) catalogue. The orbit has a 
high eccentricity and a low inclination, and is confined close 


to the Galactic plane (Dinescu et al. 19991. This cluster 


shows a split in the red giant branch. Recently, |Carretta| 
eTH!] ([20T^ analysed FLAMES GIRAFFE-fUVES spec¬ 


tra for 92 stars in the cluster and found that stars seem to 
be clustered into two discrete groups along the Na-O anti¬ 
correlation. jCarretta et al.j ( |2013 1 did not find a significant 
spread in [Fe/H], with the star-to-star variation being ~ 
0.1 dex and the mean [Fe/H] ~ -1.2. Kayser et al.|(|2008 ( 
found a clear bi-modality in CN in NGC 362 and [Worley] 
& Cottrell ( 2010| found homoge neity in s- and r-process 
abundances. Marln-Franch et al. (20091 used deep and ho¬ 


mogeneous photometry from HST to derive an age of 10.3 
Gyr for NGC 362. We selected high probability members of 
NGC 362 from the RAVE catalogue by taking into account 
the stars’ RVs, proper motions and location within the tidal 
radius of the cluster. In the next section we explore the main 
stellar parameters for these objects. 


4..3.1 Abundances, temperatures and gravities 

Most of the stars identihed as the highest likelihood mem¬ 
bers for this cluster have measured metallicities [Fe/H] ~ 
-1.7. We also hnd stars with [Fe/H] = -1.2 and -0.7 at the 
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NGC5139 



0.0 
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- 1.0 
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Figure 15. Distance vs. [Fe/H] for the members of NGC 5139 {ui Cen). In the top panel the distances derived in |Binney et al.] | |2014[ | 
are shown; in the bottom panel the distances are from [Zwitter et al.| | |2010t . Blue dots are stars with a large proper motion relative to 
uj Cen; triangles indicate stars affected by crowding. The vertical red line indicates the distance of the cluster from [Harris 1< [T9^ (2010 
edition). 
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NGC 362 



210 220 230 240 
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Figure 16. Left panel: Radial velocity-abundance plot for RV candidate members which fall within the tidal radius of NGC 362. The 
dashed lines indicate the nominal RV and [Fe/H] for the cluster reported in the literature. The blue dots are stars with large proper 
motions with respect to the measured proper motion for the cluster. Right panel: The metallicity distribution function for all RV candidate 
members which fall within the tidal radius (in red) and for those which also have proper motion measurements consistent with that of 
the cluster (black line). 
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nominal RV for NGC 362 (see Fig. 161. Blue dots in Fig. 16 
are stars with RVs similar to the cluster but they have large 
proper motions with respect to NGC 362. From the [Fe/H] 
distribution function in the right panel of Fig. most of 
the stars have [Fe/H] between -1.5 and -1.8. As noted above, 
high-resolution spectroscopy of members of this cluster has 
not yielded a significant spread in [Fe/H], although |Pila- 
chowski 119811 derived a mean [Fe/H] = -0.9 from several gi¬ 
ants while Shetrone fc Keaiie| ( |2000[ ) obtained [Fe/H] = -1.33 
from 12 giants in this cluster. The large range in metallicity 
for candidate NGC 362 members from RAVE - in particu¬ 
lar the apparent clump at [Fe/H] ~ -1.7 - is thus somewhat 
puzzling. Barring a significant spread in metallicity not de¬ 
tected in previous work, contamination by field stars and/or 
bright stars from the SMC would seem to be the most likely 
explanation for the observed abundances in our sample. 

In Fig. we show a CMD and temperature-surface 
gravity diagram of the candidate members of NGC 362, 
overplotted with BASTI isochrones for [Fe/H] = -1.3 dex 
I Carretta et al. ||2009 1 and an age of 10.5 Gyr (e.g. Marm 
Franch et al. 20091. In the left panel, five stars lie in the 


region of the RGB; of these five, four stars also fall near the 
giant branch of the isochrones in the right panel. However, 
as with NGC 3201 and cu Gen, the position of these stars 
relative to the isochrones would suggest that surface gravi¬ 
ties in the RAVE database might be underestimated for the 
given metallicity and age of the cluster. 

Table|^lists the stellar parameters derived in the RAVE 
survey database for the NGC 362 candidates. Unfortunately, 
there are no [a/Fe] measurements in the RAVE chemical 
abundance catalogue for these objects. 


We adopt the distance modulus given in Harris (19961 


(2010 edition), (m-M)v = 14.83 using a E(B-V) = 0.05 for 
this globular cluster. Fig.[^shows a few stars clearly outside 
of the isochrones for the metallicity and age reported for 
this cluster, suggesting that these objects are halo field stars 
(or members of the SMC) with a RV similar to NGC 362. 
High resolution follow-up spectroscopic observations of these 
targets would help to settle this question. 

We found that all the candidate stars associated with 
NGC 362 from the RAVE survey following the selection 
methodology described in this work have a distance signifi¬ 
cantly smaller than the nominal distance for the cluster (see 
Fig. 1181). However, there are only 4 stars with distance esti¬ 


mates from Binney et al. (20141, while Zwitter et al. (20101 


estimated the distances for all the potential members. 


5 CONCLUSIONS 

We report the identification of potential stellar members 
of NGC 3201, NGC 5139 (tu Gen) and NGC 362 in the 
Radial Velocity Experiment-Fourth Data Release (RAVE- 


DR4) survey ( Kordopatis et ah] 20131 using primarily the 
precise radial velocity derived from the RAVE spectra, for 
which 80% of the stars in the catalogue have (trv < 5 km 
s“^. These three clusters have large systemic RVs (|RV] > 
100 km s“^), making them ideal for relatively uncontami¬ 
nated radial velocity selection from the bulk of the stellar 
disc. We also used proper motions (where available) and the 
tidal radii of the three clusters to make a reasonably robust 
selection of globular cluster membership. Once candidates 


were identified, we used them to test the precision and ac¬ 
curacy of the stellar parameters derived for the stars in the 
survey. The fact that distances and ages are known with 
relatively high precision for these clusters makes them ideal 
testbeds for these kinds of measurements in the targets ob¬ 
served for the survey. 

For NGC 3201 we found a star-to-star metallicity vari¬ 
ation from -1.25 to -1.83 dex with an internal metallicity 
scatter of 0.24 dex. This is larger than the typical abundance 
errors reported in RAVE-DR4, suggesting an intrinsic metal¬ 
licity scatter in the cluster. A significant spread in metal¬ 
licity has also been seen in work based on high-resolution 
spectroscopic observations ( [Gonzalez fc Wallerstein~||1998| 
Simmerer et al. 1|2013l; however, as discussed above, other 


authors have not found evidence for a large [Fe/H] spread 
(e.g. [Covey et al.||200^ [Carretta et al. |[2009[ ), and in fact 
an explanation for this discrepancy may have recently been 
uncovered ( Mucciarelli et al.|2015 l. Nevertheless, the abso¬ 
lute metallicity values from RAVE are in good agreement 
with the values reported in these studies. Using the RAVE 
chemical abundance catalogue (Boeche et al. ]MI1 ) we found 
a mean [a/Fe] ~ 0.22, also consistent with high-resolution 
spectroscopic results. Overall there is reasonable agreement 
(to within the errors) between the RAVE-derived tempera¬ 
tures and surface gravities and the BASTI isochrones for the 
metallicities and ages reported in the literature for this clus¬ 
ter. However, if we accept the abundances from the litera¬ 
ture, the candidate members tend to lie systematically above 
the relevant isochrones for ages that are generally found for 
globular clusters (10 - 13 Gyr), suggesting that the RAVE 
analysis may be underestimating surface gravities by less 
than 0.5 dex for these metal-poor giants. Distances to clus¬ 
ters can be determined very precisely; a typical uncertainty 
in distance for a globular cluster is ~ 6% (VandenBerg et 
al.|2013 1. Binney et al. (20141 distances found distances for 


the stars in our sample ranging from 2.7 to 3.7 kpc. If these 
objects are part of the main cluster these results would sug¬ 
gest that the distances for the giants we identify as members 
of NGG 3201 are underestimated by ~ 40%. [Zwitter et al.[ 


(20101 distances for stars in our sample range from 1.3 to 6.5 


kpc; these distances have a significant scatter with respect 


to the nominal value for the cluster (D = 4.9 kpc). Binney 


et al. (20141 found an age of 10 Gyr for all the candidate 


members, which is only slightly different to the age of 11.4 


Gyr reported by Munoz et al. (20131. 


For NGC 5139 (tu Gen) we reported a large spread in 
[Fe/H] (~ 2.0 dex), ranging from near solar values to [Fe/H] 
~ -2.2 dex for the selected candidate members. This result is 
in good agreement with work done by[Sollima et al.[([2005p. 


Villanova et al. (20071 and Johnson & Pilachowski (20101 us¬ 


ing high-resolution spectroscopy. We also confirmed several 
star formation episodes in uj Gen from the RAVE metal¬ 
licity distribution function. In addition, we found a large 
spread in [a/Fe] for a given [Fe/H] using the a-elements cal¬ 
culated in [Boeche et al. ( [2011| , with [a/Fe] ranging from - 
0.23 to -1-0.32, and a mean error of ~ 0.2 dex. Large spreads 


in [a/Fe] have also been found by other studies ( Pancino 
et al.[[^02[ Villanova et al. 20071. As with NGC 3201, we 


found reasonable agreement between the Tefi and log g val¬ 
ues derived from the RAVE spectra and the selected BASTI 
isochrones. The youngest stars in the cluster are believed to 
be around 8 Gyr (Villanova et al.|[2007l. If this age limit is 
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Table 8. NGC 362 candidates selected from the RAVE data (RV, proper motions, tidal radius) and their parameters 


ID 

Te// (K) 

[Fe/H] 

log g (cgs) 

[a/Fe] J-H 

J004905.3-733108 

4000± 96 

-1.25± 0.11 

0.50± 0.31 

- 0.55 

J004217.1-740615 

4526± 101 

-2.68± 0.10 

0.12± 0.35 

- 0.38 

J005038.4-732818 

7039± 83 

-1.58± 0.09 

2.53± 0.13 

- 

J010313.6-705037 

4250± 84 

-0.75± 0.10 

5.00± 0.17 

- 0.50 

J010314.7-705115 

4669± 98 

-1.62± 0.10 

3.76± 0.19 

- 0.32 

J010314.7-705059 

4288± 75 

-1.84± 0.09 

0.15± 0.20 

- 0.53 

J010315.1-705032 

4008± 96 

-1.74± 0.11 

2.22± 0.31 

- 0.84 

J010319.0-705051 

4250± 96 

-1.25± 0.11 

l.OOib 0.31 

- 0.61 

J010335.7-705052 

4000± 96 

-1.50± 0.11 

O.OOib 0.31 

- 0.29 

J011655.9-690607 

4033± 173 

-1.21± 0.21 

0.59± 0.59 

- 0.57 


NGC362 6700 5000 3300 




J-H T,„(K) 

Figure 17. Left panel: 2MASS J, J-H Colour-Magnitude Diagram (CMD) centred on NGC 362 with a radius of 0.1 degrees. The large 
dots are the stars selected as being probable members of the cluster based on RV and location within the tidal radius; stars with large 
proper motions with respect to the cluster are in blue. Note that five stars are found on the RGB while 3 stars are clearly outliers, 
suggesting they are not members of the cluster. Right panel: Temperature-surface gravity diagram for the stars selected as probable 
members of NGC 362, overplotted with BASTI isochrones covering the range in [Fe/H] from -1.3 to -1.8 and [a/Fe] = -1-0.2 dex at an 
age of 10.5 Gyr. Four stars lie relatively close to the selected isochrones. Stars with large proper motions relative to the cluster are shown 
as blue dots. 


correct, the isochrones would suggest that the spectroscopic 
gravities for the most metal-poor stars are underestimated. 


as was also the case for NGC 3201. Binney et al. (20141 found 


most of the stars in our sample lie at distances between 2 and 
4 kpc with a peak at ~ 3 kpc, but with members spanning 


from 1 to 7 kpc. Zwitter et al. (20101 also found stars rang¬ 


ing from 1 and 7 kpc; using these distances we find a group 
of stars with a mean value around 5 kpc and [Fe/H] between 
-1.0 and -2.2 dex, in good agreement with the distances for 
to Cen found in the literature. 

An age spread from 8 to 13 Gyr between different stars 


et al. p007|). However, as discussed above, other studies (e.g. 
Sollima et al.|2005 Stanford et al.|2006 1 have found a much 


smaller or even negligible age dispersion for this cluster. Un¬ 
fortunately, RAVE data can only weakly constrain the ages 
of stars ( Binney et al.|20~ |. 

For NGC 362 we found that most of the stars identified 
in this cluster have [Fe/H] ~ -1.7 dex with a range from -0.7 
to -1.8 dex. However, high resolution spectroscopic studies 
have measured the metallicity for this cluster at [Fe/H] ~ 
-1.2, and have not found any spread in [Fe/H] ( Carretta et 
al.||2013 1. It is very unlikely that NGC 362 stars have such 


in the cluster has been reported in the literature (Villanova a wide range in [Fe/H], as colour-magnitude diagrams ob- 
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Figure 18. Distance vs. [Fe/H] for the members of NGC 362. In the top panel are plotted the distances derived in |Binney et al.| | |2014[ |; 
however, there are only 4 stars with this information. In the bottom panel the distances from |Zwitter et al.| ( |201^ are shown. Blue dots 
are stars with a large proper motion; triangles indicate stars affected by crowding. The vertical red line indicates the distance of the 
cluster from [Harris I l|1996[l (2010 edition). 
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tained for this system show very little spread (e.g. Bellazzini 
et al.|200T h There is no [a/Fe] information for these stars in 


the RAVE chemical abundance catalogue. The large discrep¬ 
ancies in both metallicity and metallicity spread between 
our sample of candidates and values in the literature suggest 
that our sample may be significantly contaminated with field 
stars and / or some of the RAVE-derived abundances may 
be erroneous. There is a general good agreement, within the 
errors, between potential members and the isochrones used 
for this cluster; however we also found that, again, using 
the metallicity and age reported in the literature for this 
cluster the isochrones suggest that the surface gravities are 
underestimated. Many candidates appear clearly outside the 
isochrones in the TeS - log g plane, supporting the idea that 
these are halo field stars or SMC stars with RVs similar to 
NGC 362. As with NGC 3201 and cu Cen, we also found that 
the distances derived in the two studies based on the RAVE 
survey ( Zwitter et al.|2010 Binney et al.|20r4 l are system¬ 
atically lower than the distance reported for this cluster in 
the literature. 


From this test of the stellar parameters, abundances and 
estimated distances of RAVE stars identified as potential 
globular cluster members, we draw two general conclusions. 
The first is that the derived stellar parameters and abun¬ 
dances are in good agreement with independent measure¬ 
ments based on high resolution spectroscopic studies from 
the literature, with the exception of surface gravities, which 
appear to be systematically underestimated relative to what 
would be predicted by BASTI isochrones for the nominal 
metallicities and ages of these clusters; in uj Cen, the cluster 


with the largest metallicity spread, this discrepancy appears 
to be strongest for the most metal-poor stars. 


The second general conclusion is that, assuming the 
candidate members are in fact associated with each cluster, 
and assuming that they are in the main body of each cluster 
(and hence at the same distance), the distances derived for 
these stars via two different methods are systematically low 
(e.g., [Binney et al.|2014p o r are distributed over a very wide 
range (e.g., Zwitter et al.|[20l0 l. While the identification of 
some of the potential members of NGC 362 is admittedly 
less secure, the preponderance of evidence indicates that we 
have found numerous genuine members of both NGC 3201 
and tu Cen in RAVE data, which, unless they are signifi¬ 
cantly extended along the line of sight, should be at approx¬ 
imately the same distance as their respective cluster. In both 
of these latter cases, the distances estimated by [Zwitter et| 
E] ( [2010[ ) appear to show a correlation with metallicity, in 
that the most metal-poor stars are assigned the greatest dis¬ 
tances. Consistent with our results, an independent analysis 
of RAVE DR4 data suggests that the distances to metal- 
poor red giants from Binney et al. (20141 may be system¬ 
atically underestimated (T. PifH, priv. comm.). Note that 


Binney et al. (20141 assumed stellar density profiles for the 


thick disk and halo while Zwitter et al. (20101 used a flat 
prior for stellar density vs. distance; in principle one might 
expect Binney et al. (20141 to prefer small distances, because 
the density of the region where a hypothetical nearby star 
is located would be higher and hence the overall solution for 
it would have a higher probability. In the case of [Zwitt^ 


et al. (20101, no constraint on stellar density and distances 
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could explain the larger distances scatter observed for the 
potential members. 

That these two general conclusions seem to be incom¬ 
patible is somewhat puzzling. As noted above, one would ex¬ 
pect systematic underestimation of surface gravities to lead 
to systematic overestimation of distances - as, at a given 
temperature and apparent magnitude, assuming a higher 
luminosity places a star at a greater distance - yet the op¬ 
posite appears to be the case with the highest probability 
RAVE globular cluster member stars. Moreover, the possi¬ 
ble underestimation of RAVE distances to metal-poor red 
giants in Binney et al. (20141 noted by PifH (priv. comm.) 
would appear to be directly opposite to the apparent trend 
with metallicity seen in the distances from IZwitter et al.l 

i Mo] ). 

This conundrum notwithstanding, the work presented 
here demonstrates the remarkable utility of survey stars 
which are members of globular clusters for testing the in¬ 
formation output by stellar survey pipelines. As astronomy 
reaps the rewards of an era of both current and future mas¬ 
sive spectroscopic surveys, the validation of the stellar pa¬ 
rameters derived from the resulting spectra is of fundamen¬ 
tal importance. 
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